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The idea to describe quantum systems within a hydrodynamic framework (quantum hydrodynamics, QHD) goes 
back to Madelung and Bohm. While such a description is formally exact for a single particle, more recently the 
concept has been applied to many-particle systems by Manfredi and Haas [Phys. Rev. B 64, 075316 (2001)] 
and received high popularity in parts of the quantum plasma community. Thereby, often the applicability limits 
of these equations are ignored, giving rise to unphysical predictions. Here we demonstrate that modified QHD 
equations for plasmas can be derived from Thomas-Fermi theory including gradient corrections. This puts 
QHD on firm grounds. At the same time this derivation yields a different prefactor, 7 = (D — 2/3 D), in front 
of the quantum (Bohm) potential which depends on the system dimensionality D. Our approach allows one to 
identify the limitations of QHD and to outline systematic improvements. 
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1 Introduction 

Dense plasmas containing quantum degenerate electrons are presently of increasing interest in many fields includ¬ 
ing condensed matter, astrophysics, warm dense matter, high energy density matter and laser plasmas, e.g. 
Despite of the experimental relevance of these systems, the theoretical desciription of quantum plasmas including 
long range Coulomb interaction remains far from completed, and a large variety of different theoretical concepts 
is being applied, including quantum Monte Carlo methods EHZl and density functional theory, e.g. oca and 
references therein. The situation is even more critical in the area of nonequilibrium plasma behavior. There 
are two main lines of research: The first are many-body theory approaches such as quantum kinetic theory and 
nonequilibrium Green functions methods, e.g. mm, linear response theory, e.g. ora and time-dependent 
density functional theory. The second line comprises particle-based computer simulations - including semiclas- 
sical molecular dynamics with quantum potentials, e.g. l3 HT9l!20l or various versions of quantum molecular 
dynamics, e.g. El- 

These simulations are often very time consuming, so it is tempting to develop simpler schemes that avoid the 
treatment of the momentum dependence entirely and resort to a much simpler fluid description. Fluid approaches 
have proven very successful in classical physics, and it is well known how to include inter-particle correlations, 
e.g. 112 1112311 . Therefore, similar approaches are of high interest for quantum systems as well. In fact, it has 
already been shown by Madelung l24l and Bohm lt25l that the Schrodinger equation for a single particle can 
also be transformed to a fluid-like form, see 12611271 for an overview. Extension to quantum many-particle 
systems date far back. Quantum hydrodynamic models for superfluid bosons were presented, among others, by 
Gross l28ll29l . Pitaevski l30ll and Biswas lUTI . Finally, the extension to plasmas—i.e. many-fermion systems— 
has been performed by Manfredi and Haas (MH) in Ref. Il32l . Similar equations have been derived by Kuzmenkov 
and co-workers 13311341 . Assuming factorization of the iV-particle wave function into equally occupied single¬ 
particle orbitals MH derived a simple set of fluid equations very similar to the one-particle case and to the case of 
bosons. Although these appear to be rather drastic assumptions, these quantum hydrodynamics (QHD) equations 
have been used in recent years in a very large number of works in quantum plasmas, for the computation of linear 
and nonlinear oscillations and waves. Applications include warm dense matter, dense astrophysical plasmas, 
electrons in metals and electron-hole plasmas in semiconductors. Furthermore, the QHD equations have been 
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extended to relativistic plasmas as well as to include the spin density as an independent variable, e.g. m- We 
also mention attempts to approximately include exchange and correlation effects Il36l . 

A crucial question is the range of validity of QHD. Unfortunately, many recent applications have ignored 
the limitations, giving rise to predictions that are in conflict with basic results of many-body quantum physics. 
As a consequence some controversy about this popular theory has emerged that concerned, among others, the 
prediction of an attractive proton-proton interaction in dense zero temperature quantum plasmas lElEMSl and 
the prediction of giant spin polarizations and spin-driven lasing EK3- Among more general criticisms we note 
the papers by Vranies et al. Il44li45l and Vladimirov et al. Il46l pointing out additional limitations. Thus a reliable 
answer about the quality and limitations of the QHD predictions is important for the field of quantum plasmas. 

In this paper we approach this issue from an entirely different side. We use, as the starting point, an indepen¬ 
dent well established continuum approach in quantum many-fermion physics: Thomas-Fermi theory (TF). We 
show that the QHD equations can be derived from TF with gradient corrections. While this gives strong concep- 
tional support for QHD, the resulting equations have a different explicit form than those proposed by MH: they 
differ in the prefactor of the Bohm term. This correspondence between QHD and TF allows for the derivation of 
improved approximations for the former based on the extensive experience accumulated in the latter. 

2 Summary of Quantum hydrodynamics equations 

The QHD equations derived by MH lf32ll consist of the continuity equation and the momentum density balance 

d t n + V_n(mi) = 0, (1) 

mn {d t + uV fl ) u = -riS7 rV - (p MB + n) , (2) 

where V is a single-particle potential energy including external and induced electrostatic (Hartree) contributions, 
V = V'ext + VI i • The closure of the inifinite system of fluid equations is made by choosing the many-body 
pressure p MB as the one of an ideal Fermi gas at T = 0, assuming locality in space and time (adiabatic local 
density approximation, ALDA): 

p MB (r, t) = ^n 0 Ep(n 0 ) ( 3 ) 

5 \ n 0 ) 

where Ep is the Fermi energy and ng is a reference density. This is analogous to classical hydrodynamics where 
/; MB = nk B T. The additional pressure n is a pure quantum effect resulting from the quantum kinetic energy 
(finite extension of quantum particles) and exists even for a single particle MM- It is related to the Bohm 
potential V B , 


f ,2 vU/U 

V7r = nVV B , Vb = - - n = A 2 , (4) 

Zm yjn 

where A is the (real) amplitude of the wavefunction. While this system (with p MB = 0) is exact for a single 
quantum particle, its extension to a many-particle system m is an approximation. The first is the assumption 
of a spatial averaging (coarse graining) over length scales of (at least) a few interparticle distances F which is 
typcial for any fluid approach 111211381 . Furthermore, aside from postulating Eq. ([3]), it is assumed that the system 
has no correlations. For a plasma this means that the coupling strength should be negligible, i.e. the ratio of the 
mean Coulomb energy to the kinetic energy, (e 2 /r) <C ( Eki n )• In the ground state this amounts to (e 2 /r) <C Ep 
[where Ep is the Fermi energy] or r s = f/a B -C 1 [a B is the Bohr radius], e.g. 1 121(27 1. Finally, to obtain 
Eq. (Q}. exchange effects are neglected; so the iV-particle wave function of the non-interacting system factorizes, 

^jv(ri,...rjv) = U.^L 1 A i (r i ,t)e~' lSi/h , (5) 

where A, and S L are real functions and, furthermore, all amplitudes are assumed equal lf32ll . ,4, ('r. t) = /Ur. t). 
The combination of these approximations makes it very difficult to assess the quality of the solutions of the QHD 
equations, for a recent analysis, see Refs. H27I1361I and references therein. 
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3 Thomas-Fermi theory approach to QHD 

A different approach to QHD is possible starting from a suitable energy or action functional of the plasma and 
performing a variational minimization. Such ideas have been used in many fields. For example, Bloch l47ll 
developed a hydrodynamical theory of the electron gas. Gross l28l and Pitaevski lf30ll derived QHD equations 
for a weakly interacting superfluid Bose condensate - the famous Gross-Pitaevski equation. The connection 
between a quantized field theory and hydrodynamic equations was given in Ref. f29l . A variational principle 
for the Gross-Pitaevski equation, based on a time-dependent Lagrange density, was studied in Ref. l48l . For an 
overview containing QHD-type equations for Bose systems, see Ref. Il49l . 


3.1 Variational formulation 

Starting point is the total energy of a one-component quantum system consisting of a mean field (Hartree, “H”) 
contribution, the energy in an external potential (“ext”), kinetic energy and a residual “exchange-correlation (xc)” 
term, written as a functional of density. 


E[n] = U ext [n] + Un[n\ +T[n] + U xc [n], ( 6 ) 

U ex t[n\ = J d 3 rV ext (r)n(r), Uu[n] = ^ J d 3 rVu([n],r)n{r). 

The ground state density profile no (r) is found by minimizing this functional, under the constraint of a given 
total particle number N = J d 3 r n(r), 


0 = 5 


E — /i J d 3 rn( r) 


with the variation yielding the chemical potential 


li[n 0 \ = V ext + Vu[n 0 ] + 


6T[r 


Sn 


+ 


6U* 


Sn 


(7) 


Stability of the ground state profile is associated with vanishing of the total force, F[no] = — V/i[no] = 0. 

The transition to hydrodynamics is straightword by allowing for perturbations of the density around the ground 
state profile and for currents (not necessarily small) that are driven by the selfconsistent force F[n(r,f)] = 
—V/t[n(r, t)}. Within the adiabatic LDA it is assumed that the functional form of /i, Eq. (0, remains the same: 

^ n (r,t) + V • [n(r,t)u(r,f)] = 0, 

d 

m—u(r, t) + mu(r, t) • Vu(r, t) = -V/i[n(r, t)]. (8) 

This hydrodynamic set of equations includes, selfconsistently via 0, all forces arising from the kinetic, Hartree, 
external, and exchange-correlation energy. It remains to choose appropriate expressions for the ideal part (kinetic 
energy) and the exchange-correlation term. Here we concentrate on the kinetic energy and neglect Uxc, since 
this is beyond standard QHD. 


3.2 Gradient corrections. Recovery of the Bohm potential of QHD 

At zero temperature the kinetic energy T can be decomposed, using a systematic gradient expansion around the 
homogeneous limit, 

T[n] = TtfM + Tgc [n]. (9) 

Here Ttf [n] is the Thomas-Fermi functional (local part), and Tqc ['«] contains the gradient corrections, with the 
result 

T[n\ = J d 3 r |f[n] n(r) + 7 + 0 [(Vn) 3 ]|, ( 10 ) 


Copyright line will be provided by the publisher 










4 


D. Michta et al.: Quantum hydrodynamics for plasmas -a Thomas-Fermi theory perspective 


where the ideal local energy density is t[n] = ^E F [n] ~ n 2 / 3 . The gradient term was first computed by von 
Weizsacker (501 for fermions and bosons in the ground state who obtained the expression (ITOt with 7 = 1 . While 
his result is correct for bosons, his calculation for fermions is not, as we briefly discuss in Sec. 13.31 The fermion 
problem was subsequently analyzed by Kompaneets ED and many others, see e.g. Ref. El, for an overview. 
The correct result of the gradient expansion for fermions is due to Kirzhnitz Ii53l who performed a gradient 
expansion of the one-particle Green function and obtained expression (ITOt with the coefficient 7 = 1/9. 

Let us now establish the connection between the kinetic energy functional (ITOl ) and the QHD equations 0. 
To this end we compute the functional derivative on the r.h.s. of the momentum balance 0, with the result 


ST T f M 

dn 

ST gc [re; 7 ] 
Sn 


= E F [n], 


= 7 


8m 


Vn 


- 2 


V 2 n 


= 7 V B [n]. 


( 11 ) 

( 12 ) 


Thus, we have established a direct connection between the energy functionals of Thomas-Fermi theory with 
gradient corrections and quantum hydrodynamics. Let us summarize the main conclusions following from (IT2l) . 

1. The functional derivative of the Thomas-Fermi term yields exactly the Fermi pressure p MB , Eq. (0, justify¬ 
ing the closure approximation of the QHD equations. 

2. The functional derivative of the gradient correction coincides, up to the coefficient 7 , with the Bohm poten¬ 
tial 0 of QHD. 

3. The Bohm potential, Vb, is directly connected with the von Weizsacker functional, (i.e. with Tqc [7 = 1] )■ 
Since the first order gradient expansion result 7 = 1 is exact for bosons, the QHD equations in the form 0 
with the Bohm potential 0 apply to (non-interacting) bosons, within the validity limits of a hydrodynamic 
description. In fact, these QHD equations are equivalent to the time-dependent Gross-Pitaevski equation 
and are widely used in the field of ultracold bosonic atoms, e.g. l49l . 

4. Another case where the Bohm potential is equivalent to the von Weizsacker gradient correction result 
(Tqc[i = 1]) is the one of distinguishable particles. An example are electrons that are tightly bound in 
atoms such that their wave function overlap is negligible. This allows for an efficient quantum trajectory- 
based description of chemical dynamics l54l . 

5. The only case where the von Weizsacker gradient correction result ( Tqc [7 = 1]) is ( at l eas t formally) correct 
for fermions is the case of a single particle or two non-interacting particles with opposite spin projection. In 
the former case the density n in Eq. i l I ()[ i is just the modulus squared of the wave function. In the latter case 
the density is the sum of the contributions from both particles. 

6 . The gradient expansion for fermions has been extended to systems of different dimensionality D by Holas et 
al. ff55ll and Salasnich |56j, who recovered the functional Tqc [ 7 ], but with the D-dependent value 7 = 

While for 3D systems this agrees with the result of Kirzhnitz ( 7 3D = 1/9), for ID systems one obtains, 
instead, 7 11:> = —1/3. For 2D systems the coefficient is zero meaning that the gradient expansion has to be 
extended to higher orders, as was done e.g. by Engel et al. El- 


3.3 Comparison of the approaches of von Weizsacker and Manfredi-Haas 

The exact agreement of the Bohm potential Vb in the QHD equations with the von Weizsacker gradient correc¬ 
tion (Tqc [7 = 1]) is striking and makes a more detailed comparison of the two approaches and their respective 
assumptions interesting. In his remarkable paper 11501 von Weizsacker considered a system of many noninteract¬ 
ing fermions (protons and neutrons in nuclei) and attempted to generalize the local Thomas-Fermi theory by the 
inclusion of gradient corrections. He used the occupation number representation, thereby properly assuring the 
Pauli principle, with the following ansatz for the single-particle orbitals in a small volume V: 

V’pM = l / ^ 1/2 [l + a(p)r]e - * pr , |p| < p F . (13) 
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where a(p)r <C 1 is assumed everywhere in V. The prefactor of the plane wave appeared to be the natural lowest 
order (in r) choice to obtain gradient corrections to the local Thomas-Fermi result which is usually derived from 
pure plain waves. Using this ansatz von Weizsacker computed the density gradient 

4 f PF o 

Vn=— J d 3 p (r, f)Rea(p), (14) 

from which he obtained the kinetic energy T[n; a] of the system parametrically depending on a(p). Minimization 
of the kinetic energy yields a variational problem for the momentum dependence of a, with the result [following, 
at least, for a power law ansatz Re a(p) = oq ■ p 1 )] a(p) = ao = const. Inserting this result in the kinetic energy 
straightforwardly recovers the result (ITOt with 7 = 1 . 

It is now interesting to make, again, contact with QHD. Indeed, using the result of this variational procedure 
gives, for the wave functions (IT3l > and the density gradient (IT4l) . 

■0 p (r) = V~ 1/2 [1 + a 0 r] e~^ pr = A p (r)e~% Sp , |p| < pf- (15) 

Vn = 2 n do- 

This yields space-dependent amplitudes of the individual orbitals that are all identical, A p (r) = A(r) = 
U -1 ' 2 [1 .. a 0 r], exactly as was assumed by Manfredi and Haas ff32l . cf. Eq. (|5j, in the derivation of the QHD 
equations. Wave functions of the type (fl5l > can, in principle, be realized for bosons in the ground state (all parti¬ 
cles occupy the orbital with p = 0). At the same, time this ansatz does not yield the correct gradient correction 
of Kirzhnitz (7 = 1/9). Therefore, even though the variational minimization was performed correctly, the ansatz 
turned out too simple an inadequate for a many-fermion system. 

Von Weizsacker himself was well aware of the substantial deviations of his result from the experimental 
data (50). He clearly pointed out the need to include correlation effects what led to the Coulomb corrections [as 
well as parity terms] in his famous formula for the total energy of nuclei. Since Coulomb corrections in nuclei are 
always small compared to the kinetic energy he correctly suspected that the problems lie in the Thomas-Fermi 
model itself. For applications to quantum plasmas or atoms this is true as well but, in addition, here Coulomb 
correlations play a much more prominent role, since they may easily be of the same order or even exceed the 
kinetic energy. Therefore, in the case of plasmas, use of the TF model (or of QHD) requires restriction to 
situations of weak coupling EE), cf. Sec. [2] 

4 Conclusions and outlook 

In summary, we have presented an analysis of quantum hydrodynamic theory from the stand point of Thomas- 
Fermi theory. We have shown that the QHD equations have the correct structure in terms of powers of the density 
and lowest order density gradients. At the same time, the pre-factor of the Bohm potential in the QHD is an order 
of magnitude too large, and it depends on the system dimensionality. For 3D systems the proper choice that is 
consistent with the energy of the homogeneous Fermi gas is 7 = 1/9. The commonly used value, 7 = 1 , applies 
to bosons in the condensate and is incompatible with fermions and the Pauli principle. Earlier QHD results for 
linear and nonlinear properties of fermions that are based on this value have to be revised Eo). 

The present approach is straightforwardly extended to finite temperatures (58l using the finite-temperature 
extension of TF-theory of Kirzhnitz ll52l and Perrot (59). Finally, it remains to remove two major deficiencies of 
QHD: the neglect of exchange and of correlations. Here, the recourse to Thomas-Fermi theory with an exchange- 
correlation potential Uxc included in Eq. <[ 6 )» gives a clear strategy how to proceed and to verify previous choices 
by Manfredi et al. (36). 
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